Kondo Physics in a Rare Earth Ion with well localized 4f electrons 



Jolanta Stankiewicz^'Q Marco Evangelisti, 1 Zachary Fisk, 2 Pedro Schlottmann, 3 ' 4 and Lev P. Gor'kov 4 ' 5 

1 Institute) de Ciencia de Matenales de Aragon and Departamento de Fisica de la Materia Condensada, 
CSIC-Universidad de Zaragoza, 50009-Zaragoza, Spam 
2 Department of Physics and Astronomy, University of California, Irvine, CA 92697, USA 
3 Department of Physics, Florida State University, Tallahassee, FL 32306, USA 
4 National High Magnetic Field Laboratory, FL 32310, USA 
5 L. D. Landau Institute for Theoretical Physics of the RAS, Chernogolovka 142432, Russia 

(Dated: March 8, 2013) 

■ Dilute Nd in simple cubic LaB6 shows electrical resistance and specific heat features at low 

temperature consistent with a Kondo scale of Tk < 0.3 K. Nd has a well localized 4f 3 J = 9/2 
Hund's Rule configuration which is not anticipated to be Kondo coupled to the conduction electrons 
in LaE>6. We conjecture that the unexpected Kondo effect arises via participation of ^/quadrupolar 
(— i ■ degrees of freedom of the Nd crystal field ground state quartet. 
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The Kondo effect underpins a broad range of correlated electron behavior, 0, Q including exotic superconductivity 
in 4f an d 5f intermetallics, [3] but is unexpected in materials with well-localized / levels. Spin interactions of 41' 
electrons with conduction states are normally either of the Heisenberg exchange or of the Coqblin-Schrieffer type. 
The former exchange is ferromagnetic and arises from the Coulomb interaction in conjunction with the antisymmetry 
of the wave function. It scales with the so-called de Gennes factor, (<? — l) 2 J(J + 1) (where g is the gyromagnetic ratio 
and J is the total angular momentum of the state), and manifests itself in the Curie- Weiss magnetic susceptibility, 
spin disorder scattering in the electrical resistivity, and the depression of superconducting transition temperatures, 
for isostructural rare earths with well localized 41 levels. Ce and Yb do not follow this scaling in most cases (as 
of course divalent Eu does not as well), because their exchange is of the Coqblin-Schrieffer type, [|| coming from 
the hybridization of the 4f levels with conduction electrons, leading to an antiferromagnetic exchange, the general 
condition for the development of the Kondo effect at low temperature. The ^/levels of Nd 3+ impurities in a metal are 
deep and only slightly hybridized. [1| The exchange interaction is expected to be of the ferromagnetic Heisenberg type 
and, for this reason, the Kondo-like features for Nd impurities in LaBg we report in this Letter are surprising. What 
is essentially an old and well understood phenomenon makes a new appearance in LaBgrNd which is both unexpected 
and needing explanation. 

The rare earth hexaborides crystallize in the CaBg simple cubic structure which can be pictured as a CsCl ar- 
rangement of Bg octahedra and rare earth atoms. The boron network requires two electrons to form a closed shell 
configuration, leaving one conduction electron per unit cell in the hexaborides formed by trivalent rare earths. This 
electronic description has been well confirmed by band structure calculations [alas well as extensive de Haas- van 
Alphen measurements, and agrees exactly with Hall effect measurements, [a, NdBg is a Tn = 8.6 K type 

I antiferromagnet. The 4f 3 J — 9/2 crystal field scheme for the full cubic site symmetry of the Nd ion has been 
determined by both inelastic neutron and Raman scattering [12j experiments to be a ground state T| quartet, 
with the Tg quartet and T 6 doublet lying at 135 K and 278 K above the ground state, respectively. An analysis of 
the temperature dependent specific heat using this magnetic level scheme, in addition to the lattice specific heat of 
LaBg, led to the surprising conclusion that the electronic specific heat of NdBg above Tjv equaled 90 mJ/mole-Nd 
K 2 . [13| This has prompted the investigation of the single ion properties of Nd dissolved in the isoelectronic, isostruc- 
tural host LaBg reported here. We find that Nd diluted in simple cubic LaBg exhibits electrical resistance, magnetic 
susceptibility, and specific heat features at low temperature, which are consistent with a Kondo scale of Tr- < 0.3 K. 
We argue that quadrupolar interactions between the ^/and the conduction electrons, which give rise to the formation 
of the spin singlet, are most likely responsible for this behavior. 

All measurements were performed on single crystals grown from Al flux using high purity elemental rare earth from 
the Ames Laboratory. Specific heat measurements using the relaxation method were carried out down to approximately 
0.3 K in a commercial setup for the < noH < 9 T magnetic field range. Magnetic susceptibility measurements 
were performed in a commercial superconducting quantum interference device magnetometer for temperatures down 
to 2 K, and in a homemade Hall microprobe magnetometer installed in a 3 He setup for the lower temperature range. 
The four-probe resistivity was measured to low temperature in a 3 He cryostat or in a PPMS system. Contact leads 
(25 /im gold wire) were soldered to the samples using pure indium. In our experiments, we used a dc (less than 1 
mA) or low- frequency ac current. We etched the crystals in concentrated HC1 to clean the surfaces before performing 
our measurements. 

We measured electrical resistivity, magnetization and heat capacity on Lai-^Nd^Bg single crystals for < x < 0.06. 
In order to avoid any uncertainty in the Nd content, the very same samples were used in all experiments. The Nd 
compositions in the LaBg crystals were in general quite different from the nominal charge compositions. Therefore 
the Nd concentration was determined from low-temperature susceptibility and magnetization measurements, such as 
shown in Fig.[TIa)-(b), using the crystal field scheme determined by neutron and Raman scattering in NdBg. Since the 
effective moments for the T| and Tg are quite close, nearly the same Nd composition would be determined from the 
susceptibility measurement assuming either ground state. The specific heat measurements confirmed unequivocally 
that the ground state of Nd ions in LaBg is the T| quartet as in pure NdBg. Fig. [TJc) shows that the initial 
susceptibility Xinit of our samples has no term proportional to x 2 . Neither we find such term in the saturation 
magnetization which varies linearly with the Nd composition. We conclude therefore that the concentration of 
magnetic pairs or clusters, which usually give rise to the x 2 variation of magnetization, is negligible in our samples, [Tij 
in agreement with electron microprobe results. On the other hand, the induced magnetization in the paramagnetic 
state is M = Xi-ff + X3^ 3 + where the nonlinear susceptibility \3 probes quadrupolar fluctuations as shown for 
several rare-earth compounds. [15|, [l|| The measured \3 f° r Lai-^Nd^Bg single crystals is shown in Fig. Hfd). Its 
behavior points unequivocally to the importance of the quadrupolar interactions as discussed below. 

A Kondo rise in the low-temperature resistivity p of dilute Lai-^Nd^Bg can be seen in Fig. [2j The comparison 
with the LaBg resistivity data, also shown in Fig. [51 rules out that the minimum in p at approximately 10 K arises 
from impurities in the source material. Beyond this temperature, the resistivity of diluted alloys varies with T as that 
of the pure LaBg sample. This indicates that the structural disorder is alike in all samples, and gives rise roughly 
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FIG. 1. (Color on-line) (a) - Low-temperature magnetic susceptibility for Lai-zNcLBg [x — 0.0095) single crystal. A x oc 
— ln(T) variation is shown for T < 1 K. (b) - M(H) at T = 2 K for various Lai-^Nd^Bs single crystals, (c) - Xinit/ X versus x 
at various temperatures, (d) - The measured nonlinear susceptibility X3 times T 3 at low temperatures. 
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FIG. 2. (Color on-line) Electrical resistivity for Lai-^NdzBe single crystals. The main panel shows the temperature dependence 
of the electrical resistivity for Lai-^Nd^Be single crystals in zero magnetic field. The Kondo contribution to the resistivity in 
the low-temperature region is shown in the inset. How this contribution varies with Nd content x is plotted in the inset of the 
inset. 



to the same residual resistivity. The magnitude of the low-temperature rise in resistivity Ap, plotted in the inset of 
Fig. [21 is 3 /ifkm per mole Nd. For the well studied case of Lai_a;Ce x B6 the value per Ce is 50 times larger, [17( where 
Tk = 0.86 K. [l8| The excess resistivity goes nearly as A — BT 1 / 2 which points to non-Fermi liquid-like behavior. We 
observe the rise in the resistivity for a very small Nd-doping only (x < 0.03) as shown in the inset of Fig. [2j With 
increasing x, one might expect sufficient strain in Lai-^Nd^Bg to lift the degeneracy of the ground-state quartet. If 
the four-fold degeneracy is central to the presence of the Kondo features found in these alloys, then the effect should 
strongly weaken with Nd concentration. 
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FIG. 3. (Color on-line) Specific heat for Lai-^Nd^Bg single crystals. Top panels (a) and (b) show the temperature dependencies 
of the experimental specific heat (C) for Lao.99osNdo.oo95B6 and Lao.97osNdo.o295B(i single crystals, respectively, along with C 
of the isomorphic LaB6, at various magnetic fields, as labelled. Bottom panels (c) and (d) show how the excess (magnetic) 
entropy S m per Nd ion, normalized to the gas constant R, varies with temperature for the same single-crystals as in (a) and 
(b), respectively. 



In panels (a) and (b) of Fig. [3] we show the T-dependence of the specific heat (C) of dilute Lai-a-Nd^Bg for the 
x = 0.0095 and 0.0295, respectively, in addition to isomorphic nonmagnetic LaBg. The strong dependence of C on 
H at the lowest T is brought about by the Nd ions. This magnetic contribution (C m ) superposes on the lattice 
contribution at higher temperatures, which is given by the data for LaB 6 compound. Furthermore, the electronic 
specific heat seen in LaBg (linear term C e i oc T in Fig. 3) should be expected to contribute likewise to Lai-^Nd^Bg. 
Therefore, C m (T,H) is easily obtained by subtracting C of LaBg from the total C of Lai-^Nd^Bg. Thus, by taking 
the integral / °°(C m /T)dT, the magnetic entropies (S m ) are calculated and depicted in panels (c) and (d) of Fig. [3] 
for x — 0.0095 and 0.0295, respectively. For {iqH > 1 T, C m (T,H) is well within our experimentally accessible 
temperature window. Indeed, extrapolating linearly C rn /T as T — > K yields an S m content for T < 0.3 K that is 
less than 2% of the total S m content. For H = 0, the lack of experimental C m for T < 0.3 K has been taken into 
account by matching the limiting S m at high T with the value obtained from the in-ficld data. After scaling S m to the 
proper Nd concentration, one can see that there is an entropy of R ln(4) per mole Nd involved in the excess specific 
heat of Lai-^Nd^Bg above that of LaB 6 , indicating that (i) the 4f 3 ground state is compensated by the conduction 
electrons and some other effects as T — > K and (it) the ground state is a quartet, as is Ce in Lax^xCe^Bg. [l9| 
Assuming that about half of the zero- field entropy of a Kondo system is reached at the characteristic temperature Tk, 
the development of S m in panels (c) and (d) would suggest that Tk is of order of < 0.3 K for the samples studied. 

The zero-field magnetic contributions of Nd ions to the specific heat (C m ) are shown in the top panel of Fig. [4] 
The dominant Schottky-type feature is seen to broaden and shift downward as x decreases from 0.0295 to 0.0095. By 
applying a magnetic field, these specific heat contributions gain height and get narrower, being no longer dependent 
on the Nd concentration for any field higher than I — 2 T (Fig. [3]). For /j,qH > 6 T, we clearly distinguish two separate 
anomalies in C m (bottom panel of Fig. 0}, which we associate with the two doublets into which the Tg quartet splits 
in a magnetic field. They have different g-factors, depending on the angle of the field with the crystallographic axis. 
The exact value of the (/-factors for Nd also depends on the crystalline field splitting of the ion. For Nd, and for the 
field along (100) direction, the ratio of the ^-factors is 2.94/1.92 ps 1.7, using the Lea-Leask-Wolf table for J = 9/2 
and x = —0.90. [2(| This is very close to what we obtain from a fitting of the Schottky four-level specific heat to the 
experimental LaB 6 :Nd data at = 6 and 9 T (solid lines in Fig. [J}. 

The observed behavior agrees with the following picture. Neglecting the possibility of a hybridization between the 
localized and conduction states, leaves us only with the ferromagnetic Coulomb exchange. We ascribe the Kondo-like 
phenomena to the coupling of the conduction electrons with the quadrupolar momentum of the Nd-sites. Because 
the magnetic degrees of freedom strongly couple with the quadrupolar ones, the degeneracy is lifted at T = 0. A 
sufficiently high magnetic field and increasing x break the Kondo-like coupling. 

The distance between Nd ions is not very large for the x values considered and is a random variable. Hence, in zero- 
applied field, the Nd ions are not totally independent but see each other. It is therefore expected that the interaction 
between Nd ions leads to an internal field that locally splits the T\ ground-state, according to a distribution whose 
width decreases with the increase of the Nd concentration. This interferes with the Kondo-like screening and leads 
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FIG. 4. (Color on-line)Magnetic specific heat for Lai-^Nd^Be single crystals. Top: Temperature dependencies of the experimen- 
tal zero-field magnetic specific heat (Cm) per Nd ion for Lao.ggosNdo.oogsBe and Lao.9705Ndo.o295Be single crystals, respectively. 
The same data, plotted as C m /T vs. log(T), are in the inset. Bottom: Magnetic specific heat for Lao.9705Ndo.o295Be single 
crystal between 0.3 and 20 K at various magnetic fields, as labelled. Solid lines are fits to the Schottky anomalies for the 
corresponding fields. 



to a distribution of "Kondo temperatures" : Tk differs at every site or locally in a region covering several unit cells. 
Both effects reduce the height of C m and make the "Kondo resonance" broader. Assuming that (i) for T < Tk the 
specific heat of an individual ion is proportional to T/Tk, and (ii) the distribution of Tk is reasonably flat, then 
C m /T oc ln(T /T), where T in an average of Tk (see inset of Fig. 2]). [21| Similarly, the low-temperature zero-field 
susceptibility is proportional to ln(To/T) (see Fig. da)). A distribution of Kondo temperatures can also explain such 
behavior. The susceptibility for an individual ion for T < Tk is proportional to 1/Tk- A flat distribution of Tk then 
again leads to x ^ — la(T).(2lj| This agrees with our data. The excess resistivity Ap oc T 1 / 2 over nearly a decade of 
temperature is also indicative of non-Fermi-liquid behavior. 

We should stress that the estimated value of Tjc (< 0.3 K) marks no more than a characteristic temperature scale 
(from competition between the ferromagnetic exchange and scattering on the quadrupolar momentum of Nd-ions) 
and does not define the magnetic field scale same way as it does in the s-d model. For LaB 6 :Nd, the magnetic-field 
scale is given by the field that splits the T\ quartet into two doublets. In Fig. @J the double-feature in the specific 
heat shows at poH = 6 T, which corresponds to T m 4 K. 

The low-T magnetic susceptibility data follows a Curie- Weiss law down to 2 K with an antiferromagnetic intercept 
of ~ —0.3 K. The nonlinear susceptibility sheds additional light on the origin of the unexpected Kondo feature seen 
in the Nd alloys. For noninteracting ions and within a quartet consisting of two doublets with Zeeman splitting, 
given by g\ and <?2, respectively, the linear and nonlinear susceptibilities are xi — (G?i + ffD/^^CMs/^s) an d 
X3 = — ((<?i +92 +6<7 2 <7f)/(192 x T 3 ))(fi^/kg), respectively. The free-ion X3 is always negative reflecting the negative 
curvature of the Brillouin function. A plot of X3^ 3 vs T then gives a constant and deviations from this constant are 
consequently brought about by the Kondo-like screening. An inspection of Fig. [IJd) shows that the measured X3T 3 
for LaBgiNd is both negative and, for T < 1 K, monotonically increasing with decreasing temperature. The observed 
magnitude and temperature dependence of X3 is similar for the other two samples studied. These results strongly 
favor the quadrupolar model in LaBg:Nd. 

Quadrupolar effects are strongly evident in the electronic properties of a number of rare earth intermetallics, most 
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notably in the Pr skutterudites where they appear to give rise to heavy Fermion behaviour. |22l . l23j However, one 
does not find single-ion Kondo behavior in the properties of dilute Pr dissolved in the La skutterudite analogs. [24j No 
heavy Fermion behavior has been seen in the stoichiometric Nd skutterudites. Raman experiments as well as detailed 
theoretical studies do indicate that one expects strong quadrupolar influence in the properties of NdB6 . [H, HI, HI] 
Note that for Lai -xCe ^Bp , strong Kondo features exist for all Ce concentrations and for large x antiferroquadrupolar 
order has been found. |28[ 

It is established that quadrupolar degrees of freedom play a fundamental role for Ce and Nd ions with T 8 ground- 
state; they manifest themselves in first place through interactions among the ions. Impurities strain the lattice and 
the local lattice distortions couple to the quadrupolar, but not the spin, degrees of freedom of the ^/electrons. The 
quadrupolar intersite interactions are probably the main reason for why it is hard to observe the single ion limit in both, 
LaBg:Nd and LaBg:Ce. The induced quadrupolar distortion could also lead to an interaction between the states 
and the conduction electrons. There is no consensus about the origin of the interactions. Kubo and Kuramoto [25[ 
successfully described the excitation spectrum of NdBg using nearest-neighbour intersite exchange and quadrupolar 
interactions. A different approach emphasizing crystalline fields was proposed by Uimin and Brenig. [27| For CeBg, on 
the other hand, quadrupolar interactions between sites, 
model, (3(| H3 and a detailed group theoretical study 



29] the RKKY interaction arising from the Coqblin-Schricffer 
32| have been presented. The well-localized nature of the 4f 



electrons in Nd compared to the considerable hybridization in Ce, suggests that the interaction mechanisms with the 
conduction states could be different for these two ions in the LaB 6 matrix. 
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